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General formalism for describing dynamics of modulated beams along linear beamlines is devel- 
oped. We describe modulated beams with spectral distribution function which represents Fourier 
transform of the conventional beam distribution function in the 6-dimensional phase space. The 
introduced spectral distribution function is localized in some region of the spectral domain for nearly 
monochromatic modulations and it can be characterized with a small number of typical parameters 
which we choose to be the lowest order moments of the spectral distribution. We study evolution 
of modulated beams in linear beamlines and find that characteristic spectral parameters transform 
linearly. The developed approach significantly simplifies analysis of various schemes proposed for 
seeding X-ray free electron lasers. We use this approach to study several recently proposed schemes 
and find the bandwidth of the output bunching in each case. 
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I. INTRODUCTION 

Free electron lasers (FELs) are the only currently avail- 
able sources capable of generating coherent X-ray radi- 
ation. Presently several hard X-ray FELs (XFELs) are 
operated or under construction . These light sources 
are designed to work in the Self- Amplified Spontaneous 
Emission (SASE) regime @ which amplifies shot noise 
within the FEL bandwidth. As a result, the output X- 
ray pulses have limited longitudinal coherence which may 
be undesirable for some applications. The XFEL per- 
formance can be strongly improved if it is seeded with 
coherent signal. Seeded FELs generate more powerful 
radiation with smaller bandwidth and their undulators 
are shorter. That allows to design cheaper light sources 
having significantly larger brightness. 

Growing FEL mode couples radiation with electron 
motion. Therefore, FEL can be seeded either with nar- 
rowband radiation Q or with electron beam modulated 
at the X-ray wavelength (8Hl0|. The second option 
looks attractive since electrons are charged particles and 
they interact with electromagnetic fields unlike radiation 
which weakly interacts with materials at X-ray frequen- 
cies. Large variety of available beam optics elements in- 
troduces a lot of options for transforming modulation 
and controlling parameters of the output beam. A num- 
ber of different schemes for generating modulated beams 
were recently propos ed |8rJl2| and some were verified ex- 
perimentally [La - fla ]. All these schemes utilize the same 
principle: first the beam is modulated and then the beam 
modulation is transformed in the following beamline to 
produce bunched beam. Conventional analysis of each 
scheme relies on following electron trajectory in the phase 
space. Such an analysis is somewhat complicated and 
does not provide intuitive understanding of the scheme 
performance when the beamline parameters are changed. 
Moreover, such an approach does not allow one to pre- 
dict the bandwidth of the resulting bunching if the initial 
modulation is not perfectly coherent. As a result, it is 
not clear whether the proposed schemes are suitable for 



FEL seeding. In this paper we address this drawback 
and develop general formalism quantifying modulation 
parameters and describing evolution of bunched beams 
along an arbitrary linear beamline. This approach will 
simplify the analysis and will serve as a powerful tool for 
designing FEL seeding schemes. 

The modulated beam has fine scale features in the 
phase space distribution. Low order moments of the 
beam distribution such as rms beam sizes, rms energy 
spread, etc. cannot capture these small scale variations. 
As a result, one needs to introduce high order moments 
to characterize the distribution function. That signifi- 
cantly complicates the analysis since the number of pa- 
rameters characterizing the phase space distribution in- 
creases by many orders of magnitude. Alternatively, one 
can consider the 6-dimensional (6D) Fourier transform 
of the beam distribution function which introduces spec- 
tral distribution function. Such a spectral distribution 
is well localized in some region of the spectral domain 
and it can be characterized with a small number of pa- 
rameters such as the wavevector of modulation and its 
bandwidth. These parameters capture the main proper- 
ties of the modulation similar to the rms beam quantities 
characterizing its envelope. This approach is very similar 
to description of laser pulses which are commonly char- 
acterized with the carrier frequency and bandwidth (low 
order moments of the spectral power distribution). The 
spectral distribution changes when the electron beam 
passes through beam optics elements. As it will be shown 
bellow, one can easily find how the spectral distribu- 
tion evolves along the beamline. That allows finding the 
change of the rms spectral parameters and eventually de- 
scribing the evolution of the bunched beam. 

The paper is organized as follows. In Sec. |TT] we de- 
velop general principles for describing dynamics of mod- 
ulated beams. First, we introduce spectral distribution 
function describing beam distribution in the spectral do- 
main and derive equation describing its change along lin- 
ear beamlines. Then we introduce low order moments 
of the spectral distribution and find how they evolve 
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along the beamline. Finally, we find several invariants 
relating spectral properties of modulated beams which 
remain constant under arbitrary linear symplectic trans- 
forms. In Sec. lIIII we describe the main methods for mod- 
ulating beams and find the resulting spectral distribu- 
tion function which is generated in these modulators. In 
Sec . 11 VI we apply developed formalism for analyzing com- 
mon schemes for XFEL seeding. We demonstrate that 
developed theory significantly simplifies analysis and es- 
timate the bandwidth of the output bunching in case of 
significant random phase noise of utilized laser pulses. 



II. GENERAL DESCRIPTION OF 
MODULATED BEAMS 

Any electron beam can be described as an ensemble 
of electrons occupying some phase space volume. Each 
electron can be fully described by its position in the 6D 
phase space 



C(s) = (x,p x ,y,p y ,At,-AE), 



(1) 



where x and y are the transverse coordinates of electron 
in respect to the reference trajectory, p x and p y are the 
corresponding canonical momenta, At is the deviation of 
arrival time to position s along the beamline, AE is the 
deviation of the particle energy from the average bunch 
energy. 

The entire bunch as an ensemble can be described with 
the distribution function /(£) which characterizes elec- 
tron density in the 6D phase space. We anticipate de- 
scribing modulated electron bunch which manifests as 
fine scale variations in the phase space. Describing such 
a beam with rms quantities will require finding high or- 
der correlations in order to capture short scale variations. 
Therefore, it is more practical to describe modulated 
beams with the entire distribution function rather than 
with a very large number of high order moments. 



method. Each characteristic describes electron trajectory 
in the phase space and generalizes Newton equations for 
an arbitrary choice of canonical variables £: 



as 



(3) 



In this paper we neglect that particles interact with each 
other (Hamiltonian does not depend on the distribution 
function). We also limit our analysis to linear beamline 
optics. Such a system can be described with Hamilto- 
nian H which quadratically depends on the phase space 
coordinate £. The Hamiltonian can be represented as a 
quadratic form defined by the symmetric matrix H = Ti. T 



'1 

ds 



5 s /(C,s) + (V/) T JHC = 0, 



(4) 
(5) 



Under this assumption, Newton equations ([3]) become 
linear 



dC 

ds 



JHC 



(6) 



and they can be solved using, for example, Magnus ex- 
pansion. The final and initial coordinates of each particle 
transform linearly by conventional symplectic transform 
matrix R. Then the formal solution of Vlasov equation 



C(s) = R(s,s )C(s ), 
f(C,s) = f(R- 1 (s : s K,s ) 



(7) 

(8) 



This solution represents the Liouville theorem which 
states that the particle density is conserved along the 
trajectory in the phase space. 



B. Spectral distribution function 



A. Vlasov equation 

Along the beamline electrons interact with electro- 
magnetic fields which satisfy Maxwell equations. That 
implies that the forces can be described with Hamiltonian 
-ff (C, s) which describes the entire dynamics. The evolu- 
tion of the distribution function satisfies Vlasov equation 
which can be considered as a continuity equation in the 
phase space 

f = d s f(C,s) + [f,H} = 0, (2) 
as 

where [/, H] = (V/) T J(Vif) is the Poison bracket, J 
is the unit block-diagonal antisymmetric symplectic ma- 
trix, V is the 6D gradient in the phase space, and super- 
script T stands for transposition. 

Vlasov equation is a hyperbolic partial differential 
equation and it can be solved using characteristics 



As mentioned above, the distribution function /(C) of 
the modulated beam has small scale variations and it can- 
not be characterized with a small number of typical pa- 
rameters. However, these variations are quasi-periodic. 
Therefore, the Fourier spectrum of the distribution func- 
tion is well localized and it is more uniform than the 
phase space distribution function. It can be used then 
for description of modulated beams. This approach is 
similar to methods used in Optics where laser pulses are 
conventionally described with their spectra rather than 
the actual field distribution in space. Unlike Optics, the 
beam is characterized with the distribution function in 
the 6D phase space and any component of the electron 6D 
position can change along the beamline. Therefore, one 
has to consider 6D Fourier transform of the distribution 
function which introduces spectral distribution function 

/ k (k, S ) = J f(C,s)e- lkT <d 6 C (9) 
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The spectral distribution function fully describes the 
beam since the phase space distribution can be recovered 
through inverse Fourier transform of the spectral distri- 
bution. At the same time, the spectral distribution func- 
tion reflects important characteristics of the beam which 
are not evident from the distribution function /(C) i n the 
phase space. For example, spectral distribution function 
along longitudinal axis ks = k z is equal to the beam 
bunching factor, b(k) = J f(C)e~ lk ^d 6 C = jk(k = k 5 fc). 

Description of modulated beams with the spectral dis- 
tribution function is beneficial only if its change along 
the beamline can be described with simple enough equa- 
tions. Otherwise, the model will be too complicated for 
practical use. Bellow we find the exact equation de- 
scribing the change of the spectral distribution but first 
we show that its dynamics is simple enough. Consider 
a single harmonic of the spectral distribution function 
/k(k, so) = S(k — k') which corresponds to the distri- 
bution function /(£, s ) oc exp(ik' T £). The phase of the 
modulation <j> = k T C linearly depends on the phase space 
coordinate C- Therefore, the phase of the modulation 
remains linear under any linear transform of the phase 
space C RC, which is the case of linear beamline op- 
tics. Therefore, a monochromatic modulation transforms 
into another monochromatic modulation having different 
wavevector, k' — > k". As a result, different harmonics of 
the spectral distribution do not mix with each other in a 
linear beamline and spectral distribution function is self- 
similar along the beamline. This property indicates that 
the equation describing evolution of the spectral distri- 
bution function should be simple enough. 

We find the precise equation describing dynamics of 
the spectral distribution function by taking Fourier trans- 
form of Vlasov equation ^) assuming quadratic form of 
the Hamiltonian described by Eq. (jj). Using equality 
tr(J'H) = one can end up with the following equation 
after some straightforward algebra 

^ = d s f k + [f k ,H u ] k = 
as 

= dsfk + (V k ./k) T H Jk = 0, (10) 

#k(k,s) = --^k T JH(s)Jk. (11) 

This equation has exactly the same form as the orig- 
inal Vlasov equation (|5|) and can be interpreted then 
as Vlasov equation for the spectral distribution func- 
tion. The dynamics of the spectral distribution function 
is symplectic and it is fully described with Hamiltonian 
H(k, s) defined by Eq. (fTT|). 

According to Liouville theorem the spectral distribu- 
tion function does not change along the trajectory in k 
space, and therefore, it can be considered as an ensem- 
ble of non-interacting quasi-particles. Each quasi-particle 
represents one spectral harmonic of the spectral distribu- 
tion. These quasi-particlcs do not vanish or born in time, 
they do not interact or mix with each other, and their en- 
tire dynamics manifests as motion along trajectories in k 



space. The only fundamental difference between descrip- 
tions in the phase space and spectral domains is that 
the spectral distribution function can be complex unlike 
the phase space distribution function which is real and 
positive. 

Each harmonic of the spectral distribution function re- 
mains as a single harmonic along linear beamline. This 
property indicates that linear beam optics cannot be 
used for generating high order harmonics of modulation. 
This observation does not contradict performance of the 
HG HG j| or EEHG schemes. As it will be shown in 
Sec. lIII 51 high order harmonics are generated in the mod- 
ulator in which the phase space of the beam transforms 
nonlincarly. The following chicanes do not create addi- 
tional harmonics in the spectral domain but only trans- 
form imposed modulation into longitudinal bunching. 

In order to solve spectral Vlasov equation (TIU)) we note 
that it is hyperbolic partial differential equation, same as 
conventional Vlasov equation ([5]). Therefore, it can be 
solved using characteristics method. Each spectral com- 
ponent remains constant along the following trajectory 
in k space 

dk , 
z = HJk. 12 

ds 

Characteristics equation (fT2)l for each spectral har- 
monic is similar to characteristics equation ([5]) describing 
trajectory of each individual particle. Therefore, these 
linear equations can be solved using the same techniques. 
At the same time, these two solutions are related to each 
other since corresponding equations contain the same 
matrices. One can note that two arbitrary solutions of 
Eqns. (|6|) and (IT21) do not evolve independently and their 
product is invariant along the beamline 

^=0, (13, 
ds 

Using solution for particle trajectory (|7|) and noticing 
that relation (|13[) holds for any arbitrary initial parti- 
cle position C( s o)j one finds the solution for the spectral 
characteristics 

k(s) = R- T (s,s )k(s ). (14) 

This relation describes the change of each spectral com- 
ponent wavevector. The solution of the spectral Vlasov 
eqaution (fT0|) reads then as 

/k(M) = fr(R T (s,s )k,s ). (15) 



C. Spectral averages 

The beam is fully described with its spectral distribu- 
tion function /k(k, s). At the same time, few average 
parameters of the spectral distribution function may suf- 
fice to characterize its main features. We introduce the 
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spectral averaging of some variable <?(k) as 
Jg(k)|/ k (k, S )| 2 rf 6 k 

f|/ k (k, S )|^k ' 



.9(k)( S ) = 



(16) 



Note that we use overline notation 7 to characterize spec- 
tral averages in order to distinguish them from the aver- 
ages in the phase space domain defined as < • >. The 
denominator in Eq. (j!6[) remains constant in linear beam- 
lines and serves as the normalization factor so that the 
quantity itself and its spectral average have the same di- 
mensions. 

In case of narrowband modulations the spectral distri- 
bution function is localized in some region of the spectral 
domain. Such a distribution can be well characterized 
with few low order moments. We introduce the spectral 
average wavevector of modulation k(s) and the second 
order spectral correlation matrix 



B(s) = (k-k)(k-k) T . 



(17) 



This matrix describes rms spreads of the spectrum dis- 
tribution function /k(k, s). The diagonal elements de- 
scribe the modulation bandwidths along the correspond- 
ing axes. Therefore, we call matrix B as the "bandwidth 
matrix". The bandwidth matrix defines the rms ellip- 
soid in the spectral domain similar to the rms envelope 
ellipsoid defined by the beam matrix E =< QQ T >. This 
rms bandwidth ellipsoid qualitatively shows noise distri- 
bution along different dimensions. 

In some cases the spectral distribution function 
/k(k, s) can be localized in several well separated regions 
of the spectral domain, e.g. the spectrum of the beam 
modulated with the laser pulse consists of several local- 
ized harmonics as will be shown in Sec. IIIIBI Vlasov 
equation (fTOj) for the spectral distribution function states 
that different trajectories in the spectral domain do not 
cross which does not allow well separated spectral do- 
mains mixing with each other. In this case the spectral 
average quantities can be introduced for each localized 
domain individually. 

Using solution (fT5|) for the spectral distribution func- 
tion along the beamline, one can easily find transforms 
of the spectral averages 



k(a) = 
B(s) 



ir T ( s , So )k( So ), 

= R- T (a,a )B(a )R-\a,a ). 



(18) 
(19) 



Note that the wavenumber of modulation transforms 
linearly along the beamline. This property indicates that 
scaling the wavenumber of modulation without changing 
the wavevector orientation results in the same scaling 
of the output modulation wavenumber. This property 
can be used, for example, in the masked dogleg [16| or 
Emittance EXchanger (EEX) [l4[ schemes to generate 
bunches with different spatial spacing. This task can 
be accomplished by using several masks with different 
spacing between the slits. At the same time, the optics 
required to transform imposed modulation into longitu- 
dinal bunching remains the same for different masks. 



The bandwidth matrix B(s) transforms exactly in the 
same way as the inverse beam matrix E _1 , 

E(s) =< CC T >= R(s, .s )E(so)i? T (s, s ). (20) 

This property becomes evident if one considers unmodu- 
lated beam with some envelope. In this case the average 
modulation wavenumber is close to zero and the band- 
width matrix describes the rms wavevector spread of the 
spectral distribution which is inverse proportional to the 
rms beam sizes. Moreover, if one considers the beam with 
6D Gaussian distribution, / oc exp(— ^ T E _1 C/2), then 
the bandwidth matrix is related to the rms beam ma- 
trix as B = E _1 /2 (this property can be proved through 
direct calculation of integrals in Eqs. © and (jTB")) by 
transforming integration variables according to Eq. (|A1[) 
which diagonalizes matrices). This relation also holds 
if arbitrary monochromatic modulation is imposed on 
the beam with 6D Gaussian envelope since this modu- 
lation changes the average wavevector but does not af- 
fect the bandwidth matrix. At the same time, if the 
imposed modulation is not monochromatic and has some 
bandwidth due to noise, then the resulting bandwidth 
matrix deviates from Y,~ 1 /2. This example shows that 
describing modulated beam with the bandwidth matrix 
is useful only in the presence of significant noise when 
|| J — 2JBE|| 3> ||/||. Otherwise, the modulation can be 
considered monochromatic and it is more convenient to 
describe the beam with its phase space envelope and a 
single wavevector of modulation. Introducing the band- 
width matrix for the monochromatic modulation is not 
feasible since it is strongly related to the beam matrix E 
and does not provide additional information. 

Note that the bandwidth matrix is a useful concept for 
describing the beam quality even if the beam is not modu- 
lated. While the rms beam matrix E describes large scale 
features of the beam such as its rms sizes, energy, and an- 
gular spreads, the bandwidth matrix B describes small- 
scale fluctuations in the distribution function. There- 
fore, the bandwidth matrix describes the beam homo- 
geneity and it is more homogeneous at smaller values of 
||7-2BE||. 

The modulation wavevector and the bunch envelope 
transform independently from each other since their 
transforms depend only on the transform matrix R but 
not on each other. Identical envelopes transform in the 
same wave regardless what kind of modulation is imposed 
on the bunch. Identical modulations imposed on bunches 
with different envelopes also transform in the same way. 
This property significantly simplifies the task of design- 
ing beam buncher since it can be designed independently 
from the rest of the beamline. This beamline section 
should solve the problem of transforming imposed mod- 
ulations into required modulations and this optics does 
not depend on rms beam quantities such as its sizes and 
cmittances. Once designed, this block can be attached to 
any accelerator and it will perform equally in any regime 
of the machine. 
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D. Modulation invariants 

The bandwidth matrix is a positively defined symmet- 
ric matrix which transforms along the beamline through 
symplectic matrix R T JR = J since it represents Hamil- 
tonian dynamics. This transform is similar to the trans- 
form (|20|) of the beam rms matrix E. Therefore, the 
invariants which hold for the beam matrix (l7j are also 
applicable for the bandwidth matrix 

det(.B) = inv, (21) 
tr(BJ) 2 " = inv, n = 1,2, ...dim(B)/2. (22) 

These invariants can be used to introduce parameters 
similar to eigen-emittances for the beam envelope [i~7l |. 

Other than that, one can note that the change of the 
spectral averages is fully described with the transform 
matrix R. This is not surprising since the transform ma- 
trix defines the mapping of each particle in the phase 
space, and therefore, describes the entire beam dynam- 
ics. Therefore, spectral averages do not transform fully 
independently from the beam rms quantities. Combin- 
ing Eqs. (fTS)) and (|20|) one can easily find the following 
scalar quantity which remains constant along arbitrary 
linear beamline 

k T (s)E(s)k(s) = inv. (23) 

This quantity can be interpreted as the overall phase of 
the modulation across the bunch envelope. To illustrate 
this invariant one may consider imposing longitudinal 
bunching on the beam and its further compression. In 
this case, the number of modulation periods within the 
rms beam length remains constant during compression, 
i.e. the overall phase of modulation is invariant. 

Additional invariants can be found by taking into ac- 
count that transform matrix R describes Hamiltonian 
dynamics and, therefore, is symplectic. This condition 
limits possible changes of the beam matrix and modula- 
tion wavenumber. Using Eqs. (fT5|) . (j2U)) . and taking into 
account that transform matrix R is symplectic one can 
construct the following invariants 

k T ./(S" 1 J) 2,l+1 k = k T ( JST 1 ) 2 '^ 1 Jk = inv, (24) 
n = 1,2, ...dim(£)/2. 

There is infinite number of invariants but the number of 
functionally independent ones is equal to the dimension 
of space describing beam dynamics, i.e. it is equal to 
3 in general case. Moreover, these invariants are not 
functionally independent from Eq. ([23]). 

Invariants (|24|) can be understood using eigen- 
cmittance concept as described in Appendix |Al They 
state that the overall phase of modulation across each 
eigen- phase plane is preserved. This property can guide 
attempts for beam bunching using transverse masking. 
It is of particular practical interest to consider beams 
which eigen-emittances are recovered as x-, y-, and z- 



cmittances since this scenario corresponds to the bright- 
est beam with the fixed phase space volume. Transverse 
masking modulates beam only in a single eigen- phase 
plane (which coincides with x- phase plane) in case of un- 
corrected beam. According to invariants (|24|) the mod- 
ulation imposed in one eigen- phase plane remains in the 
same eigen- phase plane along the beamline. This indi- 
cates that transform of x-modulation into z-modulation 
should be accompanied with transform of x- eigen- phase 
plane into z- eigen- phase plane. Therefore, one should 
use EEX optics [lU to achieve this goal. Otherwise, 
the modulation can be smeared out at significantly large 
emittances as observed in Ref. [lj| . 

The bandwidth matrix transforms exactly in the same 
way as the inverse rms beam matrix. Therefore, the in- 
variants involving the bandwidth matrix are similar to 
invariants (|23|) and ((24)) found for the beam matrix 

kV 1 k = mu, (25) 
k T J(B J) 2ri+1 k = inv, n=l,2,3. (26) 

Invariant described with Eq. (|25p can be interpreted as 
the conservation of relative bandwidth Sk/k along the 
beamline. In case of homogeneous modulations, invariant 
([25]) is exactly the same as invariant (|23|) since the band- 
width matrix is proportional to the inverse beam matrix 
and the relative modulation bandwidth is roughly equal 
to the inverse number of modulation periods per bunch. 



III. BEAM MODULATORS 

In this section we describe ways to impose modula- 
tions on the beam. Currently, there are two known ways 
for modulating beams, namely transverse masking and 
imposing spatially dependant energy modulation. 



A. Transverse masking 

The beam can be modified by passing it through the 
mask. The mask serves as a filter which either absorbs 
some particles or strongly changes their phase space co- 
ordinates so that they can be filtered downstream the 
beamline. Significant change of the phase space coordi- 
nates for some group of particles while keeping intact the 
rest of the beam requires presence of non-Hamiltonian 
forces such as incoherent damping, scattering, ionization, 
etc. The simplest mask can be envisioned as a perfo- 
rated foil installed transversally to the bunch propaga- 
tion. Electrons passing through the foil strongly scatter 
and are lost downstream while electrons passing through 
the opening remain intact. These masks modify the 
bunch transversally and they are widely used in accel- 
erator technology. Modifying the distribution function 
along other phase space coordinates would require dif- 
ferentiation of particles based on their position in phase 
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space: different ionization cross-section at different ener- 
gies, triggering interaction externally at a certain time, 
etc. 

Without addressing particular physical mechanism of 
masking, we consider that the bunch distribution func- 
tion can be changed by the mask as 

/(C) = AC)M(C), (27) 

where / is the distribution function of the bunch before 
the mask and M is the the mask imprint. The spectral 
distribution function of the resulting beam /k(k) can be 
presented as a convolution of the initial and the mask 
spectral distribution functions /£(k) and Mk(k), respec- 
tively 

/ k (k) = f° * M k = J /^(k')M k (k - k')d 6 k', (28) 

This equation provides a semi-qualitative description 
for the modulated beam spectral properties. We consider 
both the beam and the mask spectra to be well-localized. 
If their bandwidths can be neglected i.e. their spectra 
can be well approximated with (5-functions, the wavevec- 
tor of the resulting modulation can be approximated as 

k[/k]«k[/£]+k[M k ]. (29) 

Note that using perforated foil as a mask modulates the 
beam in the transverse plane. For example, a set of wires 
aligned along y and placed periodically along x would 
shift the average wavevector of modulation along k x . 

Taking into account finite bandwidths of the initial 
beam and the mask results in the finite bandwidth of the 
resulted modulation. In some cases (Gaussian or uncor- 
related spectra) the spectral bandwidths add in quadra- 
ture to each other. In general case, the bandwidth of the 
resulting beam cannot be related to the bandwidths of 
the initial beam and the mask without knowledge about 
the spectral distribution functions. It can be just noted 
that if one of the bandwidth is much larger then another, 
then the bandwidth of the convolution is approximately 
equal to the largest one. 

B. Energy modulation with laser pulse 

In this scheme the beam is modulated by applying spa- 
tially dependent energy modulation Q. Such a modula- 
tion can be imposed by passing a beam through the undu- 
lator and its interaction with resonant radiation. Inter- 
action of the laser pulse with electron beam results in the 
change of the particle energy E — > E — SE smtp(z). The 
resulting beam distribution function f(z,E) relates to 
the distribution function of unperturbed beam f°(z,E) 
as 

f(z,E) = f(z,E + 6Esm l p(z)), (30) 

where for simplicity of description we omitted transverse 
variables. 



Ideally, the laser is considered to be monochromatic, 
i.e. f(z) = kz. In order to extend our analysis we 
consider that the laser has finite bandwidth caused by 
random fluctuations in its phase (and consider the laser 
amplitude to be spatially uniform). As a result, the laser 
phase has linearly growing and random components. In 
order to characterize such a laser, we consider average 
quantities of its spectrum 

MQ= feW-^dz, (31) 



k v = k[(p k ], = ( fc - k v) 2 i<Pk]. (32) 

The spectral distribution function of the modulated 
beam described with Eq. ([30)) can be related to the spec- 
tral distribution functions of the initial beam and the 
spectral distribution of the laser phase ipi,. After some 
straightforward algebra one can find that the beam spec- 
trum consists of infinite number of harmonics 

oo 

h{k z ,k E )= E fL n) (k z ,k E )= (33) 

n— — oo 

oo 

= J n (\kESE\)f^(k z ,kE) *<Ph{k z ) * — * <fk(k z ), 

V v ' 

n— — oo 

n 

where J n is the n-th order Bessel function. The schemat- 
ics of the spectrum is presented in Fig. [TJ 

The form of the spectral distribution function de- 
scribed with Eq. (|33j) determines the main parameters of 
the modulation spectrum. The spectrum consists of the 
infinite number of harmonics. The spectral distribution 
function along longitudinal wavenumber k z is described 
with n times convolution of the laser spectral distribu- 
tion ipk(kz)- As discussed in Sec. A[) the spectral 
properties of such a distribution cannot be found in gen- 
eral case. However, if one assumes that the laser phase 
noise is stationary and Gaussian (the laser spectrum is 5- 
correlated which can be assumed if the laser bandwidth 
is mainly defined by the noise rather than its envelope 
and the spectrum statistics is Gaussian) then the spec- 
tral properties of the beam modulation can be found as 

M/k°] = M/k]+^ (34) 
[/ k " } ] = B*„*. [/£] + \n\SkJ, (35) 

i.e. the carrier frequency scales linearly with the har- 
monic number and the bandwidths add in quadrature. At 
the same time, if the laser phase is not random then the 
modulation bandwidth increases much faster with har- 
monic number. For example, if one considers chirped 
laser pulse then the modulation bandwidth would scale 
linearly with harmonic number, Bk,,k z = Bk t ,k t [fk] + 
n 2 Jfc^ . This effect can be proved by expressing n— th or- 
der convolution as (pk(k z )* ■■■* ( Pk(k z ) = J e lnv ^~' lk '' z dz, 
which indicates linear scaling of the chirp with harmonic 
number. 
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Unlike the longitudinal spectral distribution along k z , 
the energy spectral distribution along k E is equal to the 
product of the initial beam spectrum and the modulation 
spectrum which is equal to the Bessel functions. The 
product of these two spectra does not vanish as long as 
both distributions are significant in the same domain. 
The schematics of the beam spectrum resulted from its 
energy modulation with laser is shown in Fig. [TJ Dashed 
area schematically shows the energy bandwidth caused 
by the initial beam envelope. The harmonics having sig- 
nificant amplitude within that bandwidth are present in 
the spectrum and high order harmonics are heavily sup- 
pressed. 



where erg is the rms energy spread of the original beam 
envelope. Note that this condition is harder to satisfy 
for modulations at higher harmonics which agrees with 
scalings of HGHG scheme Q . 

The energy bandwidth of the resulting distribution 
{Bk E ,k E element of the bandwidth matrix) is not well 
defined since the spectral energy distribution consists 
of several spikes which are described with oscillating 
Bessel function. One cannot consider them individu- 
ally since they are not well separated in the spectral 
domain. However, one can characterize the spectral dis- 
tribution bandwidth with its curvature at the maximum, 
B kE .k E 0-5fa(k z ,k E )/dl E kE fk(k Z7 k E ). For low or- 
der harmonics, /x„i < 5E/a E , one finds 




B 



(n) 



fi n) (k z ,k B ) 



1 

1 0.618n 2/3 



Mni ~ n 2 {5Ef 



{5Ef 



n <Ll ( 37 ) 



Eq. (|33|) shows that the non-diagonal terms of the re- 
sulting bandwidth matrix Bk zl k E = Bk E .k, are equal to 
zero if these terms were zero for the original beam enve- 
lope before the modulation was applied. This approxi- 
mation is typically valid for beam without energy chirp. 
Then the bandwidth matrix of the modulated beam can 
be approximated as 



B[f (n) 



Bk^kAfi] 



\n\8kl 











(38) 



FIG. 1: (Color online) Spectrum of modulation produced by 
spatial variation of beam energy. Dashed area schematically 
shows the energy bandwidth of the beam envelope. 

In order to find the energy wavenumber {kg ) max 
which maximizes each harmonic one needs to maximize 
the spectral distribution function 



^-(J n (\k B SE\)f°(k z ,k E )) 

Uk E 



0. 



(36) 



This equation can be solved assuming certain beam en- 
ergy distribution, otherwise the problem cannot be solved 
analytically. Typically the beam is assumed to have 
Gaussian energy distribution (so as its energy spectrum) 
[1, H, [f9j]. The main features can be recovered consid- 
ering beam energy bandwidth to be very large, i.e. the 
beam energy spread is considered to be small. Then the 
maximum of n-th harmonic is not affected by the energy 
bandwidth of the original beam envelope and it is located 
at 

max w n n \/8E, where [i n \ is the first maximum 
of n-th order Bessel function, J' n (^ n i) = 0. Therefore, 
the modulation is not heavily suppressed by the finite 
beam energy spread if the modulation is significantly 

large (k E ) max < {B kE .k E [f^}) 1/2 or SE > \i nl o E ~ na E , 



Note that this expression describes characteristic spec- 
tral bandwidth rather than exact spectral average value. 
Therefore, transforms and invariants for the bandwidth 
matrix found in Sec. Ill CI and III Dl should be used with 
caution. 



IV. TRANSFORMS OF BEAM MODULATION 
IN VARIOUS SCHEMES 



In this Section we illustrate the developed formalism 
describing evolution of modulated beams in the spectral 
domain. We consider common schemes for XFEL seeding 
which rely on beam modulation by laser and transform- 
ing imposed modulation into longitudinal bunching at 
small wavelengths. As described in Sec. IIII Bl the beam 
dynamics occurs only in the longitudinal phase space, 
C = {t,-AE) T . Therefore, we limit our analysis to the 
2D distribution function in the longitudinal phase space 
and its spectrum. 

Conventional schemes for creating beam bunching 
from the laser-induced energy modulation rely on two 
types of beam optics elements, namely chicanes and RF 
cavities introducing energy chirp. These elements trans- 
form the wavevector of modulation in 2D spectral domain 
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as k = R T k°, where the transform matrices for these 
elements have the following form: 



1 



R 



cavity 



l 

a 1 



(39) 



As a result, the wavevector of modulation transforms by 
these elements as 



k z 

k E 

k z 
k E 



= R 



chi< 



= R 



-T 

cavity 



k° 
K E 



k° z 

k° 
h E 



(40) 



(41) 



It is easy to follow these transforms in the 2D spectral 
plane similar to what was presented in Fig. [TJ Any given 
scheme can be represented with a diagram which shows 
trajectory of the modulation wavenumber between the 
initial and final values. The following notations on these 
diagrams will be used: 

• Consider the laser induced modulation imposed 
on the beam. This modulation has spectral dis- 
tribution function which peaks at the wavevector 
k = (2im/\, (kg ) max ) T , where A is the laser wave- 
length and (tg') mal can be found from Eq. (|36|) . 
We will focus on the evolution of the vicinity of 
this spectral component since it has the largest har- 
monic current and we will mark this wavevector 
with a cross. 

• Each beam optics element transforms the modula- 
tion wavevector according to Eqs. (|40f and ([41~j). 
The modulation wavevector changes its k E com- 
ponent when the beam passes through some chi- 
cane and it changes k z component when the beam 
passes through some RF cavity. We will illustrate 
these transforms with arrows starting from the ini- 
tial and ending at the final wavevectors of mod- 
ulations. Note that chicanes and cavities will be 
represented as arrows parallel to the spectral space 
axes. Also note that chicanes are typically consid- 
ered to have positive dispersion £ > 0. Therefore, 
they will be represented with downward arrows in 
the right-hand side of the diagram (positive har- 
monics numbers, k z > 0) and upward arrows in the 
left-had side of the diagram (negative harmonics 
numbers, k z < 0). 

• Some schemes require imposing laser-induced mod- 
ulation several times along the beamlinc. As dis- 
cussed in Sec. IIII Bl the modulation does not vanish 
only in the vicinity of k E where both the beam en- 
velope and the modulation have significant spectral 
distributions. In these schemes the harmonic cur- 
rent imposed by the first laser and transformed in 
the beamline can serve as an envelope in the follow- 
ing modulator. We will mark this mechanism by a 
dashed arrow parallel to k z axis (modulation by the 



second laser does not vanish only in the vicinity 
of k E of already existing modulations). Not that 
this arrow refers to the nonlinear modification of 
the beam spectrum unlike solid arrows which refer 
to the linear transforms in which different spectral 
harmonics do not mix with each other. 

All the schemes considered bellow are designed to 
transform initial modulation into longitudinal bunching. 
Therefore, the linear beamline optics is aimed to trans- 
form the wavevector of the imposed modulation into 
k = (fc 2 ,0) T , i.e. to eliminate any energy modulation. 



A. High Gain Harmonic Generation (HGHG) 

In this scheme the energy modulation imposed by the 
laser pulse in a short undulator section is linearly trans- 
formed by the following chicane Q. The schematics for 
the transform of the modulation wavevector is presented 
in Fig. [5] Considering the final modulation to be longitu- 
dinal bunching (i.e. k E = 0), one can easily recover the 
chicane strength which maximizes the output bunching 
factor 



(k [n) \ 

X^E /max 

27m/A 



A 



2ttSE' 



for 



5E 

OE 



> n > 1 



(42) 



As discussed in Sec. IIII Bl and illustrated in Fig. [TJ 
the modulation at high harmonics is heavily suppressed 
because the (k E )max ~ n/SE component of the im- 
posed modulation wavenumber increases with harmonic 
number and eventually does not fit within the energy 
bandwidth of the beam envelope. That results in re- 
duced amplitude of the spectral distribution as follows 
from Eq. (|33|) . The largest harmonic number which 
is present in the beam spectrum can be estimated as 
6E/a E . 



k E >\ 



X 



FIG. 2: (Color online) Schematics for the modulation trans- 
form in the HGHG scheme. 

The bandwidth matrix can be found using Eq. (|38[) . 
Assuming that the bandwidth is mainly determined by 
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the laser phase noise, one obtains the following band- 
width matrix of the modulation at the end of the HGHG 
scheme 



B 



(n) 

HGHG 



- ^chTcane^t/k ]^ 



nSk? 



hicane 



|(43) 



The bandwidth matrix transformed by the chicane is not 
diagonal anymore. This indicates that the rms band- 
width ellipse is not aligned along k z and k E axes. There- 
fore, the extent of this rms bandwidth ellipse along k z 
axis (which defines the bunching bandwidth) can be sig- 
nificantly smaller than the Bk zt k z bandwidth matrix el- 
ement. The bandwidth of the resulting bunching can 
be found in the limit of large harmonic number, n 3> 
1,/ini ps n + 0.81™ 1 / 3 



(») 



HGHG 



1 



HGHG 



{nk v ,Qf 
-1 



= I§(i + 3 .23n 4 / 3 §] . (44) 
ktp \ k^p J 

This estimate indicates that the relative bunching band- 
width at n— th harmonic is roughly y/n times smaller 
than the relative bandwidth of the laser with J-correlated 
phase noise. The second term in Eq. (|4"4")) shows ad- 
ditional noise filtering due to the presence of the chi- 
cane. However, this term is close to unity for low or- 
der harmonics (n ~ 5) and relatively narrowband lasers 

4\ 



Sk-^/kp 



B. Echo Enabled Harmonic Generation (EEHG) 



As discussed in Sees. IIII Bl and IIV A[ the main lim- 
itation on generating modulation at high harmonics of 
the laser wavelength comes from the limited overlap be- 
tween the energy bandwidths of the beam envelope and 
the laser induced modulation. As illustrated in Fig. [TJ 
the modulation at high harmonics is located at high val- 
ues of (fc^ ) max ~ n/SE which eventually fall outside 
the energy bandwidth of the beam envelope. As a result, 
the amplitude of the modulation rapidly drops as indi- 
cated by Eq. (|3"3"|) . This drawback can be compensated by 
increasing the amplitude of the imposed energy modula- 
tion. Such an approach increases modulation amplitude 
but also results in the increased rms beam energy spread 
which may reduce FEL performance. Alternatively, the 
beam envelope can be modified in such a way that its 
energy bandwidth is shifted to the domain in which high 
order harmonics are produced. This is precisely the sce- 
nario which takes place in the EEHG scheme |{|. In 
this scheme the beam is first modulated with the laser 
pulse. The following strong chicane is used to increase 



the k E component of the nx = — 1 harmonic of modula- 
tion. Then the second modulator is used and nx = — 1 
harmonic of the already imposed modulation serves as 
the beam envelope for the secondary modulation. As 
a result, harmonics are generated in the domain of in- 
creased kE values, i.e. high order harmonics with n<i 3> 1 
can be generated. The schematics of the EEHG scheme 
is presented in Fig. [3] 




FIG. 3: (Color online) Schematics for the modulation trans- 
form in the EEHG scheme. 



Condition for the overlap of two modulations along kE 
coordinate determines the first chicane strength (for the 
optimized performance of the EEHG scheme at a given 
harmonic number). The strength of the second chicane 
can be determined from the condition that the final mod- 
ulation should be recovered as the longitudinal bunching, 
i.e. ks = after the second chicane. Taking these con- 
ditions simultaneously results exactly in the same equa- 
tions as presented in Ref. @. Note that the first chicane 
is used to increase kE component of the beam modula- 
tion. This process results in strong energy modulation 
of the beam, i.e. creating energy bands within the beam 
envelope which are observed in simulations Q. 

The resulting bunching is linearly transformed from 
the modulation produced in the second modulator. The 
energy bandwidth Bk E ,k E is mostly determined by the 
energy bandwidth of rii = — 1 harmonic generated in the 
first modulator for n = n\ + 712 ^> 1, 5E\ ~ 8E2, which 
follows from Eq. (|37[) . At the same time, the spatial 
bandwidth of the resulting modulation B^ z ^ z increases 
every time the laser modulation is applied. Finding the 
precise bandwidth matrix of the resulting modulation is 
complicated since the bandwidth matrix is not diagonal 
before the second modulator due to presence of the first 
chicane. However, in the typical regime of small laser 

bandwidth, Sk^/k^ — > 0, the bandwidth matrix remains 
mostly diagonal throughout the beamline and one can 
estimate the bunching bandwidth generated in EEHG 
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scheme as 

(«) 



.2 



5k 



2 5k% 



EEHG 



n 2 IT 2 ' 

• "CD 



8E X 
~5E~ 2 



5k 2 

=£ « I- (45) 

krn 



This scaling indicates that the relative bunching band- 
width reduces as at large harmonic numbers, same 
as for HGHG scheme. 



C. Compressed Harmonic Generation (CHG) 

CHG scheme was proposed to generate short wave- 
length bunching through longitudinal bunch compression 
|1G | . This effect can be understood by considering the in- 
variant (|23[) . If the beamline is designed to compress the 
bunch, the bunching wavelength decreases to keep the 
number of modulation periods within the bunch fixed. 
Such a scenario can be realized by passing the beam 
through the RF cavity which introduces energy chirp and 
the following bunch compression in the chicane. As it was 
noticed in Ref . [h| , such a design requires strong RF cav- 
ities to introduce significant compression without smear- 
ing out the resulting bunching factor. Alternatively, the 
resulting bunching smearing can be avoided by decreas- 
ing the energy modulation amplitude but then the re- 
sulting bunching factor rapidly drops down same as in 
HGHG scheme. This effect can be understood by con- 
sidering modulation transform in the spectral domain. 
RF cavity changes the longitudinal component of the 
beam modulation wavevector k z = k® — ak® E according 
to Eq. PT|). The change of the longitudinal wavenum- 
ber is proportional to kg component of the modulation. 
That implies that the largest change of the modulation 
longitudinal wavenumber can be reached at small laser- 
induced energy modulations since fc^ oc 1/SE. At the 
same time, the amplitude of the laser-induced energy 
modulation cannot be much smaller than the rms beam 
energy spread to provide large amplitude of modulations 
as follows from Eq. (|33| . Therefore, the largest change 
in the longitudinal wavenumber of strong enough modu- 
lation is limited to Ak z ~ a/as, which requires strong 
RF cavities for significant compression of modulations 
(increasing the longitudinal wavenumber by an order of 
magnitude). 

To overcome this limit of the CHG scheme, the sec- 
ond chicane was introduced in the beamline |10j . The 
first chicane increases kg component of the modulation 
so that the following cavity would result in larger shift of 
the longitudinal wavenumber k z as follows from Eq. (|41l) . 
Once the longitudinal wavenumber of modulation is in- 
creased, the second chicane is used to recover modulation 
as purely longitudinal bunching. The second RF cavity 
can be optionally added to the beamline to eliminate the 
energy chirp from the beam. This cavity does not change 
the longitudinal wavenumber of modulation (kg = 
for longitudinal bunching) and therefore, does not affect 
beam bunching once it is generated. The schematics of 



the described scheme is presented in Fig. [4] The sequence 
of the RF cavities and chicanes duplicates the design de- 
scribed in Ref. [l(| • The final bunching is recovered from 
the n = — 1 harmonic of the laser modulation in order to 
utilize chicanes with positive energy dispersion. 




FIG. 4: (Color online) Schematics for the modulation trans- 
form in the CHG scheme. 



The scheme showed in Fig. |4] is designed to trans- 
form initial modulation having a given wavenumber k° = 
{—kip, (4"')max) T into longitudinal bunching having the 
wavevector k = (Mfc v ,0) T , where M is the modulation 
compression factor. Such a transform can be accom- 
plished with a linear beamline consisting of two chicanes 
and the RF cavity accelerating the beam at X-point 



ch i,cane2-**-cavity 



chi< 



!r T k° 



(46) 



Two equations relating components of the modulation 
wavevector provide relationships between three param- 
eters of the beamline (dispersions of two chicanes and 
the RF cavity strength). Assuming that the amplitude 
of energy modulation significantly exceeds the rms beam 
= fj.il/SE, one finds 



energy, i.e. (k^) 



6 



6 



M + l _ /zn 
—ai SE 
M+l 1 



M 



-o l 



(47) 
(48) 



Note that the required cavity strength is negative to 
ensure positive chicane dispersion. Also note that the 
strengths of both chicanes are inversely proportional to 
the amplitude of the RF field. Therefore, this scheme 
is a subject for parameters study to ensure that all the 
elements can be made linear enough. 

The bunching bandwidth can be estimated using in- 
variant ()25l) . The imposed modulation is mostly longi- 

tudinal, k v a 2 3> (J, 2 (6 E) 2 so the compression can 
also be considered as mostly longitudinal. Then the rela- 
tive bandwidth of the modulation is preserved along the 
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linear bcamlinc and one can estimate CHG bandwidth as 
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(49) 



The relative bunching bandwidth mirrors the relative 
bandwidth of the laser pulse having phase noise and it 
does not depend on the beam compression factor M. 
Note that this scheme requires smaller laser bandwidth 
compared to HGHG and EEHG schemes in order to 
achieve the same bandwidth of the resulting bunching. 



V. CONCLUSIONS 

In this paper we developed formalism for quantitative 
description of modulated beams dynamics. This formal- 
ism is based on introduction of the spectral distribution 
function which is equal to the Fourier transform of the 
beam distribution function in the 6D phase space. The 
introduced spectral distribution is localized in some re- 
gion of the spectral domain which makes possible to char- 
acterize it with a limited number of parameters such as 
the average wavevector of modulation and the bandwidth 
matrix. The key advantage of the proposed formalism is 
based on the observation that any given spectral compo- 
nent of the distribution remains as a single component 
under linear transforms indicating that the topology of 
the spectral distribution remains the same along linear 
bcamlinc. As a result, the characteristic parameters of 
the spectral distribution transform linearly. 

The developed formalism has two main advantages 
compared to alternative approaches based on trailing 
coordinates of each individual electron in the phase 
space. First, the formalism quantifies the main parame- 
ters of modulated beams which introduces suitable met- 
rics characterizing modulated beams. Different compet- 
ing schemes can be easily compared to each other using 
this metrics. Second, the developed formalism describes 
the change of the spectral parameters along arbitrary lin- 
ear beamline. The spectral parameters transform linearly 
and their change is solely described with the conventional 
beam transform matrix. As a result, the developed de- 
scription of modulated beams dynamics is full and self- 
consistent. 

This description estimates the bandwidth of the result- 
ing modulation which allows to determine whether some 
proposed scheme satisfies requirements for a given appli- 
cation. As an example, we considered several laser based 
schemes proposed for XFEL seeding and found the out- 
put bunching bandwidth if the laser has random phase 
noise. 

Note that description of the beam in the spectral do- 
main of the phase space is a useful concept for wide class 
of problems in beam physics. It can be applied for solv- 
ing homogeneous linear problems (20| or analyzing linear 
instabilities. 
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Appendix A: Interpretation of modulation invariants 



Invariants (|24[) of the modulated beam dynamics can 
be interpreted using the eigen-emittance concept 17 1. 
According to Williamson's theorem [2lj . any real sym- 



metric positively defined matrix E can be diagonalized 
with a symplectic matrix A, i.e. A T JA = J, through 
the transform 



A(s)H(s)A T (s) = A, 



A 2n - 



l,2n 



-1 ' A2„ : 2ri — £ 



hij = for i ± j, (Al) 
(A2) 



(n) 



In case of the rms beam matrix E =< ££ T >, 
are called the "eigen-emittances" and they remain con- 
stant along the linear beamline. Williamson's theorem 
can be interpreted in two ways. The first interpretation 
states that any arbitrary beam matrix can be diagonal- 
ized throughout appropriate beam optics which trans- 
form matrix is equal to R = A. Alternatively, matrix A 
can be considered as canonical transformation of canon- 
ical variables C ci „ = A£ in which the beam matrix be- 
comes diagonal. Williamson's Theorem states that the 
transform matrix A can be chosen in such a way that the 
diagonal terms come in pairs, i.e. A2 n -i,2n— l = A-2n,2n- 
However, that requires renormalization of the canonical 
conjugate variables so that they have equal dimensions. 
As a result, the elements of the beam matrix lack of clear 
physical interpretation such as rms bunch sizes or angu- 
lar spread. Alternatively, we chose transform matrix A 
which only decouples the eigen- phase planes but does 
not necessarily diagonalize the beam matrices associated 
with these planes 



A(s)E(s)A T (s) = diag (s (n) (s) 
det fs^(s) 



» 



mv, 



(A3) 
(A4) 



where E< n )(s) is 2 x 2 beam matrices corresponding to 
n-th eigen- phase plane. 

The concept of eigen- phase planes provides simple pic- 
ture for bunch dynamics since particle motion in different 
eigen- phase planes is decoupled from each other. Note 
that decomposition into the eigen- phase space coordi- 
nates is not unique and so the beam dynamics within 
a given eigen- phase plane depends on the choice the of 
eigen-coordinates. For example, if one chooses the eigen- 
coordinates £ eig = A(s)£ so that A(s) = A(sq)R(s, so)" 1 
then the beam is stationary in this rotating frame, i.e. 
<9 s E(")(s) = 0. This choice of rotating phase space may 
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be not the simplest description. For example, quad focus- 
ing and diffraction are easier to describe in a stationary 
frame. 

The eigen- phase plane concept also simplifies descrip- 
tion of modulated beams. The change of the phase space 
coordinates from the lab frame to the eigen-coordinatcs 
C cig = A(s)£ defines the wavevector of modulation in 
this frame k C i g (s) = A~ T (s)k(s). One can naturally de- 
fine the projections of the wavevector onto the eigen- 

phase planes as k^ = ((fc cig ) 2 „-i, (fc cig ) 2n ) T . These 
projections of the wavevector evolve independently from 
each other since particle dynamics within different eigen- 
phase planes is not coupled. Moreover, if one chooses 
rotating frame such as A(s) = A(so)R(s, so) -1 then the 
wavevector of modulation in this frame does not change 
along the beamline. As a result, the entire distribution 
function remains constant in this frame since each indi- 
vidual Fourier harmonic does not change. The choice of 
rotating frame which does not eliminate particle dynam- 
ics within each eigen- phase plane, (9 s E(")(s) 7^ 0, also 
implies evolution of the eigen-wavevector, d s \C' n ' (s) 7^ 0. 
However, in this case one can find the invariants 

(k^ (s)) T (s)kW ( s ) = inv. (A5) 



These invariants can be interpreted as the rms phase of 
the modulation across each eigen- phase plane area. 

The same invariants can be formulated in the lab frame 
without introducing cigen-coordinates. In this case the 
invariants read as 

k T (s) J (lT 1 (sy) 2 ™ +1 k(s) = inv, (A6) 
n = 1,2, ...dim(£)/2. 

There is infinite number of invariants but the number 
of functionally independent ones is equal to the dimen- 
sion of space describing beam dynamics, i.e. is equal 
to 3 in general case. Moreover, preservation of the rms 
modulation phase across each eigen- phase plane leads 
to preservation of the rms modulation phase across the 
entire phase space volume 

dim(E)/2 T 

(W n) ) E^kW = k T Sk = inv. (A7) 

n=l 
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